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Abstract 



This paper deals with a polymeric matrix composite material. The matrix behaviour is described 
by the modified Rabotnov's nonlinear viscoelastic model assuming the material is nonlinear viscoelas- 
tic. The parameters of creep and stress-relaxation kernels of the model are determined. From the 
experimental data related to kernels approximated by spline functions and by means of the method of 
, weighted residual, the formulas for the determination of viscoelastic parameters are derived. 



;> ; 1 Introduction 

Most of polymeric matrix composite materials are characterized by a nonlinear viscoelastic behaviour, even 
at moderate loading levels. A long term behaviour modeling these materials requires the determination 
of viscoelastic characteristics. For the hereditary- type model the methods for determining viscoelastic 
characteristics play a role of great importance and reduce to the establishment and identification of creep 
and stress-relaxation kernels. The past few years, many authors studied the nonlinear viscoelastic matrix 
behaviour modelled as per Schapery's constitutive law. See [T| and references therein. Approaches based on 
coupling Schapcry model and viscoplastic model proposed by Zapas and Crissman are used to predict a long 
term material behaviour by homogenization [2]. In the phenomenon approach of Goldcnblat-Kopnov [3], 
the hereditary model was developed for building long term strength phenomenon tests for anisotropic 
composite materials. This nonlinear model is based on the coupling of plastic potential of Mises-Hill for 
anisotropic materials and Iliuchin motion approach. In this case, taking into account the peculiarities of 
material mechanical properties, the viscoelastic parameters can be determined on the basis of long term 
strength tests of material. 

In a series of papers [IHS], it was proposed the identification methods for Schapery model using 
the uniaxial tests of creep-recovery. This approach presents two disadvantages: first, it is too difficult to 
reproduce exactly the theoretical creep recovery tests, and secondly, one generally deals with the creep 
recovery tests separately. 

In this paper, in order to avoid these difficulties, we propose a method of hereditary kernel identi- 
fication for nonlinear viscoelastic materials as follows. The experimental data of the kernels are provided 
by achieving independant tests whose number is equal to the multiple of the number of kernel unknowns. 
The precision of parameter values of the viscoelastic model requires an efficient method for determining 
hereditary parameters of nonlinear viscoelastic materials whose matrix behaviour is modelled as per modi- 
fied Rabotnov's hereditary integral equations with creep and stress-relaxation kernels, taking into account 
the experimental parameter data for small time values. The experimental data of kernels are approxi- 
mate by spline functions, and, by means of the method of weighted residual, we obtain formulas for the 
determination of viscoelastic parameters. 



1 



2 Mathematical model 



For nonlinear viscoelastic materials, the matrix behaviour is modelled by the modified Rabotnov's nonlinear 
hereditary-type integral equation [TU] 

<Po(e(f)) = <t) + A / K (t - r)a(T)dT (1) 
Jo 

a(t) = Me(t)) ~ X I R (t - T)M<r))dT (2) 
Jo 

where K and R represent the hereditary creep and the stress-relaxation kernels defined, respectively, by 

( T) -(t-r)«^ r[(l-a)(l + n)] (3) 

A is a parameter; a and /3 stand for the kernel parameters, t is the time and T the function defined by 

poo 

T(z)= / e- x x z ~ x dx. (5) 



e and a are, respectively, the deformation and the stress-relaxation functions depending on t. The scalar 
function <po, depending nonlinearily on e, is well approximated by a power function as follows: 

Me) = ^ (6) 
The coefficients H and q arc provided by uniaxial traction experimental data. 



3 Hereditary kernel identification 

The choice of efficient method for the kernel parameters determination depends essentially on the physical 
context. So, with respect to stress-relaxation kernel R of model (|1|2|) . from the condition e(t) = e(0) = 
const, the time derivative of equation ([2]) gives 

_ 1 q da(t) 

Since the above expression of the kernel R is proportional to the stress-relaxation velocity da J^ with the 
coefficient of proportionality remaining constant for s(t) = const, we can determine the kernel parameters 
of the model described by ([TJ and @ by the derivation of the experimental curves of the kernel R. 
With respect to the creep kernel K of this model ([I} and (0) , the physics imposes to satisfy the condition 
a(t) = er(0) = const that provides the time derivative of equation ([T]) as follows: 

1 Hjeit)}"- 1 de(t) 
K(t)= X a dT (8) 

As this expression of the kernel K is proportional to the creep velocity ^jp-, with the coefficient of 
proportionality remaining variable for a(t) — const, the above approach is no longer applicable, i. e. we 
cannot determine the kernel parameters by derivation of the experimental curves of kernel. In this case, the 
kernel parameters are provided by using the similarity condition between the stress-relaxation isochronic 
curves and the creep curve of the model described by ([T]) and ([2]) [TO] : 

W) (e,0) = (l + G(t)M(e,t) (9) 



2 



where 1 + G(t) is the similarity function defined as 

1 + G(t) = 1 + A f K(r)dr. (f 0) 

Jo 

Here, ipo(.) is the function represented by the deformation curve and ft(-)-, the function represented by 
the creep isochronic curves at every instant t. 

Setting S{t) = 1 + G(t), the condition © can be differentiated with respect to t, yielding the 
deformation velocity as follows: 

ds(t) f d^[S(t)a] \dS(t) 

— ; — = c i \ i — t — ; — * cr = const (11) 

dt \ d[S(t)a] J dt y ' 

which, taking into account the relation ([5]). allows to write 

4 Creep data approximation and method of weighted residual 

In this section, we use spline functions to approximate the creep kernel experimental data. Then, by 
method of weighted residual, we obtain formulas for the viscoelastic parameter value computation. 

4.1 Approximation by spline functions 

The similarity of curves representing functions (ft(-) and tpo(-) is considered in the plane (e, (p) for each 
fixed deformation level Ei,i = and for the parameter tj,j — 1, • • • ,n. Then data ipo(ei,0) = —ef 

can be approximated by the quantity (p t (e^ ,tj). For this purpose we consider the functional II defined as 

n(^-W) = [M^-Sjit^tie^)] 2 , (13) 

i=l 

where Sj (t) is the similarity mean function. Then the functional minimum is reached for 

Sj(t) = E -f o(£t ' 0) ^ (£ f j) . (14) 

In order to obtain the best approximation of the discrete data of the similarity mean function, wc choose 
the spline functions defined as [TT] 

Sj(t) = ~Aj + Bj(t - tj) + - t 3 f + Dj(t - tjf (15) 

where Aj, Bj,Cj, Dj are the coefficients, j = 1, ...,n. By substituting the equation (fT5|) into equation 
(H2J), the experimental data K(tj) can be estimated by the functions Kj(t,q) defined on [tj,tj + i] as 



Kj(t, q) = Bj + 2Cj(t - tj) + 3Dj(t - i,) 2 , j = !,-■■ ,n (16) 



with the coefficients 



Bj = Kit,), (17) 
hj-i[Mj - tij-i) 

3D, = K ^)- K ^) Dl = , (i 9 ) 

hj = tj+x-tj. (20) 
The functions Kj(t, q) obtained by approximation are presented in Table 1. 



3 



3 




1 ) 

b j 


2Cj 




Kj{t,q) 


1 





3750 








OTC A 

3750 




5 


3500 


1 AA 

— 100 


1 A 

— 10 


OdAA 1AA/-t r\ 1A/-t E \ 2 

3500 — 100(£ — 5) — 10(£ — 5) 


o 
O 




dzoO 


1 A A 

— i4y 


1 A A O 


onrrA A A(\{ ± '7^ 1A /) O / -A T \ 2 

JzoO — 149(£ — rj — 10, 4z(£ — 7) 


4 


1 n 
10 


zyOO 


— lo7 


— D.OO 


OA A A 1 7/ -* 1A\ Ci OP. ! 4- 1 A \ 2 

zyOO — lo((t — 10 j — o.oo(£ — 10) 


5 




2600 


— 167 


—6.82 


Or*AA AP^ll-t- 1 n\ /? oo/-t 10\2 

zoOO — lo7(£ — 12) — b.o2(t — 12) 


6 


15 


OO C A 

zzoO 


1 OA 

— 130 


A OO 

— 4.dz 


OO C A 10A/-t i r\ A oo f X 1 r" \ 2 

zzoO — lo0(£ — 15) — A.62(t — 15) 


7 


17 


1 AAA 

1900 


— 186 


—5.47 


1 AAA 1 Of/j- A CT A T ( 4- 1 T\2 

1900 — 186(r — 17) — 5.47(r — 17) 


o 



on 
30 


1500 


on 
—39.3 


— U.65 


1 cnn on o/^ on\ n at^(+ oc\\'A 
1500 — 39.3(1 — 30J — 0.65(1 — 30) 


n 

y 


7n 
/U 


115U 


in or: 


n no 
— U.U9 


115U — lz.ZO(t — ' Uj — u.uy(t — /UJ 


10 


80 


900 


-27 


-0.17 


900 - 27(i - 80) - 0.17(t - 80)* 


11 


100 


750 


-8.3 


-0.04 


750 - 8.30 - 100) - 0.040 - 100)* 


12 


150 


500 


-6 


-0.02 


500 - 6(i - 150) - 0.02(t - 150)* 


13 


250 


300 


-2.5 


-0.005 


300 - 2.50 - 250) - 0.0050 - 250)* 


14 


350 


250 


-0.6 


-0.0008 


250 - 0.60 - 350) - 0.00080 - 350)* 


15 


750 


150 


-0.34 


-0.0002 


150 - 0.34(> - 75 °) - 0.0002(i - 750)* 


16 


1050 


100 


-0.2 


-0.0001 


100 - 0.2(i - 1050) - 0.00010 - 1050)* 



Tabic 1: Approximated function Kj(t,q) 
4.2 Method of weighted residual 

The test of best approximation of data {K(tj), j = 1, . . . , n} by the functions {Kj(t, q), j = 1, 
remains the method of weighted residual seeking to minimize the residual (error) as 



n[X,q] = \ Wj (t)[K(tj) ~ **;(*,?)]} 

7 = 1 ^ > 



. ,n} 
(21) 



for a finite set of weighting functions Wj, j = 1, . . . , n defined as 



Wj (t) = { 1 
satisfying the following conditions 



K^-XK^q) 



K(U)-XK 3 (U,q) 



K(t j ) = XK j (t,q), Wj (t) = l. 



(22) 



(23) 
(24) 



n is the number of discrct values of creep kernel for t from to £* = 1050 and m, the order of difference 
moments, m — 2, 3, 4, • • ■ . 

The parameter A and the creep kernel parameter q of the equation (|22|) arc defined in two steps: 



• At the first stage we determine experimentally the initial values Ao and qo of parameters A and q in 
the equation (|22|) ; 



At the second stage, the values Ao et qo are used, yielding 

Kit^-XoKj^qo) 



K{U)-X Q K 3 {U,q ) 



(25) 



Then the weighting functions Wj can be provided on the basis of the minimal value S m i n of residual 



6 = J2{ti j (t)lK(t j )-X K j (t,q )}\ 



(26) 



where the sum decreases as m increases. 
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Furthermore, the minimum of residual 



is reached for 



n , 

A -> fi[A, <z] = E ^(i) [K(t 3 ) - XK 3 (t, q)] 



E?=i <z) 
The very small given value 7 of residual with given j 



yields 



|^)-A ^(i,go)r 



(27) 



(28) 



(29) 



(30) 



Taking into account the equations (|28jl and (|30[) , we obtain the formula for the parameter A value compu- 
tation 



A 



E 



gfa)gj(t,g) 

j= - |iffe)-A ^(t,go)| 2 



E 



j|Jf(t J )-AoJif i (t l gd) 



For instance, when t = tj, the formula pip reduces to 



A = 



S 1(1 



jKi-Ao)^)! 5 



y a 2 fa) 



= 1. 



(31) 
(32) 

(33) 



For each fixed deformation level i = 1, • ■ ■ , I and for the parameter tj, j = 1, ■ ■ • , n, the equation ([T]) 
reduces to 



a [1 + A / A',- (r, g)dr] e? - W = 0, 



H 



{l + A[-Bjtj - Cjt| + Djt?]}, a = const 



which admits the solution q e]0, q[, where q must satisfy the inequality 

4 < m- 

Therefore we obtain the formula for the parameter cf value computation as follows: 

e? — rjjq = 0, i = !■■■, I, j — 1, ■ ■ • , 16. 
Finally, we arrive at the following formulas for the computation of the parameters A and q: 



A 



E \K 



^ IK^-XoK^qo)^ 



E \K 



3=1 



t G [tj,tj+i],j = !,■•■ ,15, 



0, i = i... ,i, i = i, 



16. 



(34) 
(35) 

(36) 
(37) 

(38) 

(39) 
(40) 
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5 Concluding remarks 



The most important element for the kernel identification method of nonlinear viscoelastic models reveals to 
be the experimental data {K(tj), j = 1, . . . , n} related to the kernels which we applied the approximation 
test to. Among all known standard tests, the cubic spline method gave the best approximation and 
using the method of weighted residual to minimize the residual error, we obtained the formulas needed to 
computate the parameters A and q. 

Thus, this study allowed us to identify the creep and stress-relaxation kernels of the nonlinear 
viscoelastic materials. Using the spline functions we approximated the experimental data related to the 
kernel that permitted, with the help of the method of weighted residual, to obtain formulas useful for the 
parameters computation of the nonlinear model. 
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